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Abstract 

To investigate genetic variability and differentiation in Asian elephant (Elephas maximus) a total of 
53 individuals from Sri Lanka, southern India, northeastern India, northern and southern Myanmar 
(Burma), northern and eastern Thailand, and Vietnam were screened for sequence variation at the mi¬ 
tochondrial cytochrome b gene. An about 480 bp part of the gene was PCR amplified from hair bulbs, 
directly cycle-sequenced, and 335 bp were analysed on a direct blot sequencing device. Eight haplo- 
types, differing from one another by one to six third position transitions at a total of seven polymorphic 
sites, were found. They differed from published sequences of African elephant (Loxodonta africana) by 
at least 18, and from those of the Siberian woolly mammoth (Mammiithus primigenius) by at least eight 
nucleotide substitutions. Haplotypes separated out by nucleotide divergence into two major groups, 
which showed differences in their geographic distribution. However, there was no indication of a major 
separation between elephants from Sri Lanka and the Asian mainland, previously assumed to represent 
two different subspecies. Populations studied showed remarkable differences as to estimates of haplo- 
type and nucleotide diversity. Phylogeographic patterns of haplotypes suggest that Asian elephants 
have formed a coherent population after their Pliocene invasion of southern and southeastern Asia. 
Taking into account records on population histories, the present local differentiation can be explained 
by varying human impact on population structure and size. 


Introduction 

Asian elephant (Elephas maximus) once was broadly distributed from Persia through the 
Indian subcontinent into southern and southeastern Asia, and into China up to the 
Yangtse Kiang. However, as a consequence of habitat destruction, hunting, poaching and 
capturing by man, this species has experienced a rapid and dramatical decline in popula¬ 
tion numbers (Sukumar 1990; Kurt et al. 1995). Even though the grand total of elephants 
in Asia still ranges from 34,000 to 56,000 (Sukumar 1990), given the subdivision of popu¬ 
lations into several large matrilinear family clans and the fragmentation of habitats, esti¬ 
mates of the respective genetically effective population size (Ne) do not exceed about 
1,000 in the best case (Sri Lanka, see McKay 1973). 

Genetic management of both wild and captive populations of Asian elephant requires 
knowledge on levels of genetic diversity within and among populations, and on systematic 
relationships among populations living in different parts of southern and southeastern 
Asia. Both demands are currently only scarcely fullfilled. Population genetic data avail- 
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able to date are based on protein electrophoresis, carried out in a small number of indivi¬ 
duals, and give a quite contradictory picture as to levels of polymorphism and heterozyg¬ 
osity (Drysdale and Florkiewicz 1989; Nozawa and Shotake 1990; Hartl et al. 1995). 
Also the systematics of Asian elephant is far from being settled. Three subspecies - 
E. maximiis maximus on Sri Lanka, E. m. indicus in India and Indochina, and 
E. m. sumatranus on Sumatra - are presumed to exist (Haynes 1991), but electrophoretic 
evidence as to a major separation between elephants from Sri Lanka and from India 
(elephants from Sumatra were not studied) proved to be inconsistent (Nozawa and Sho¬ 
take 1990; Hartl et al. 1995). 

The present study aims at characterizing genetic diversity in Asian elephant as well as 
differentiation among populations from Sri Lanka, India, Myanmar (Burma), Thailand, 
and Vietnam on the basis of cytochrome b (Cyt b) sequences of mitochondrial DNA, 
PCR-amplified from hair bulbs. 


Material and methods 

Sample collection 

Hair samples were collected from a total of 53 Asian elephants living in southern Asian elephant camps 
or European zoos. Specimens originated from southern Sri Lanka (n = 2), eastern/northeastern Sri Lan¬ 
ka (n = 12), southern India (Kerala, n = 9), northeastern India (n = 3), northern Myanmar (n = 8), 
southern Myanmar (n = 3), northern Thailand (n = 5), eastern Thailand (n = 6), and Vietnam (n = 5) 
(Fig. 1). To achieve larger sample sizes, five study units were formed: 1. Sri Lanka with a continuous ele¬ 
phant population; 2. Southern India, isolated by the Arabian Sea towards the south and by extensive 
forest clearings towards the north; 3. Northeastern India and northern Myanmar, where there is a large, 
more or less continuous population in Assam and Nagaland on the Indian side, and Myitkyina and 
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Upper Chidwan on the Myanmar side; 4. Southern Myanmar and northern Thailand, since a consider¬ 
able movement of elephants occurs between those areas; 5. Eastern Thailand and Vietnam, with a con¬ 
tinuous population until recent times (cf. Sukumar 1992). 


Amplification and sequencing of the mitochondrial cytochrome b gene 

An about 480 bp part of the mitochondrial cytochrome b gene was amplified via the Polymerase Chain 
Reaction (PCR) using the primers L 14724 and H 15149 (Irwin et al. 1991). Assay conditions were 
2.5 U Taq-polymerase (Perkin Elmer), 0.2 pM of each primer, and 50 pM of each nucleotide ATP, CTP, 
GTP, and TTP in a reaction volume of 75 pi. The bulb end of a single hair was set directly into the reac¬ 
tion solution, without any prior treatment. After an initial denaturation step at 97 °C for 15 min, 39 cy¬ 
cles were carried out with denaturation at 93 °C for 1:30 min, annealing at 53 °C for 1:15 min, and 
extension at 72 °C for 1:30 min, followed by a final extension at 72 °C for 3 min. 

The PCR-products were purified using a commercial kit (Qiagen) and cycle-sequenced using the Se- 
quiTherm Cycle Sequencing Kit (Epicentre) and the Digoxigenin-labelled oligonucleotide 5'-DIG- 
GCTTGATATGAAAAACCATCGTTG-3' for 30 cycles, each with a denaturation at 95 °C for 30 s, an 
annealing at 52.7 °C for 30 s, and an extension at 70 °C for 1 min. Samples were run on a direct blot se¬ 
quencing device (Richterich et al. 1989). After blotting onto a nylon membrane, sequences were de¬ 
tected using an anti-Digoxigenin/Alkaline phosphatase conjugate (Boehringer) and the chemi¬ 
luminescent substrate CDP-Star (Tropix), following manufacturer’s instructions. 


Data analysis 

Based on 335 bp scored sequence of the cytochrome b gene, mitochondrial haplotypes were defined and 
compared in terms of pairwise nucleotide divergence (Nei and Jin 1989), using the Neighbor-Joining 
method (Saitou and Nei 1987) and the PHYLIP 3.5c computer package (Felsenstein 1993). Published 
sequences of the African elephant {Loxodonta africana, Irwin et al. 1991), and of the Siberian woolly 
mammoth {Mammiithus primigenius, Hagelberg et al. 1994) served as outgroups. A median graph of 
the relationships between mitochondrial haplotypes was determined according to Bandelt (1992). 
Mean nucleotide diversities within and between geographic subsets of the analysed samples were calcu¬ 
lated (cf. Quinn and White 1987). Standard errors of these values were estimated by bootstrap resam¬ 
pling of nucleotides for 1000 times. Haplotype diversities at different localities were calculated 
according to Nei and Tajima (1981). 


Results 

Compared with the corresponding Cyt b sequence in the African elephant, the sequence 
of the Asian elephant differed by 18 to 25 nucleotide substitutions. Excluding the first 
95 nucleotides not scored by Hagelberg et al. (1994), Asian elephant differed from both 
the African elephant and the Siberian woolly mammoth by 8 to 15 nucleotide substitu¬ 
tions. All but one substitutions were transitions. The only transversion (A <=^ C at position 
129) found separated the African elephant from the other two species (Fig. 2). In the 
53 Asian elephants analysed, a total of eight mitochondrial haplotypes were found, differ¬ 
ing from one another by one to six nucleotides at altogether seven polymorphic sites 
(Fig. 3). All seven mutations were transitions at the third codon position, without an ef¬ 
fect on the amino acid composition of the protein. 

Using the Neighbor-Joining procedure, haplotypes separated out by nucleotide diver¬ 
gence into two major clusters A and B, consisting of MAX I, II, III, V, VIII, and of MAX 
IV, VI, VII, respectively (Fig. 4). Each cluster contained one of the two most common mi¬ 
tochondrial haplotypes, MAX V (found in 18 specimens) and MAX VI (found in 12 speci¬ 
mens), respectively. Within each cluster, all other haplotypes can be derived by one to 
two transitions from these frequent haplotypes (Fig. 5). While the frequent haplotypes 
MAX V and MAX VI were distributed over the whole study area, some of the other hap- 
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lotypes (e.g. MAX I, MAX IV, MAX VIII) were found only in particular study units 
(Tab. 1). 



Fig. 2. Nucleotide sequences of part of the Cyt b gene in Asian elephants {E. m.) in comparison to pub¬ 
lished sequences from one African elephant (L. a., Irwin et al. 1991), and two Siberian woolly mam¬ 
moths (M.p.l, M.p.2, Hagelberg et al. 1994). Additional capital letters indicate polymorphic sites in 

Asian elephant. 
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Fig. 3. Sequence polymorphisms in 
the different haplotypes of Asian 
elephants. Position numbers are 
indicated by three vertical digits. 



Fig. 4. Sequence divergence among haplotypes of Asian ele¬ 
phant (Neighbor-joining tree, based on nucleotide divergence 
according to Nei and Jin 1989). The African elephant 
(Loxodonta) was used as an outgroup (according to our results 
the transition/transversion ratio was estimated to be 20/1). A 
and B refer to clusters of similar haplotypes as defined in 
Figs. 2 and 3. 



Fig. 5. Median graph according to Bandelt (1992), showing relationships among haplotypes in terms of 
nucleotide substitutions (small bars). Supposed transient haplotypes not found in this study are indi¬ 
cated by circles (see Figs. 2 and 3 for definition of haplotypes). 


Nucleotide diversity and haplotype diversity, characterizing genetic variability within 
study units, are given in table 1. Haplotype diversity was highest in Sri Lanka and in 
northeastern India/northern Myanmar. The latter area also showed the highest nucleotide 
diversity. Genetic differentiation among study units in terms of nucleotide diversity is 
given in table 2. 
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Table 1. Geographic distribution of cytochrome b haplotypes (MAX I-MAX VIII) found in Asian ele¬ 
phant. Letters in parentheses refer to major groups of haplotypes as defined in Fig. 4. Total T = total 
number of haplotypes found in the respective samples. SSL = southern Sri Lanka, NE SL = northern 
and eastern Sri Lanka, S In (NE In) = southern (northeastern) India, N My (S My) = northern (south¬ 
ern) Myanmar, N Tha (S Tha) = northern (southern) Thailand, Viet = Vietnam. Numbers in parenth¬ 
eses refer to study units created by pooling of single samples, h = haplotype diversity, n = nucleotide 
diversity (in per cent) for the respective study units (standard error in parentheses). 


Type 

S SL NE SL 

Sin 

NE In N My 

S My N Tha 

E Tha 

Viet 

MAX 

(1) (1) 

(2) 

(3) 

(3) 

(4) (4) 

(5) 

(5) 


n = 2 n = 12 

n = 9 

n =; 

3 n = 8 

n = 3 n = 5 

n = 6 

n = 5 

1(A) 

2 







II (A) 

3 



4 




III (A) 

1 




3 

1 

1 

IV (B) 

1 







V(A) 

1 4 

5 

1 

1 

3 

2 

1 

VI (B) 

1 1 

2 

1 

1 


3 

3 

VII (B) 



1 

2 

2 



VIII (A) 


2 






total T: 

2 6 

3 

3 

4 

1 2 

3 

3 

h 

0.84 

0.67 


0.80 

0.75 

0.65 


7Z 

0.91 

0.77 


1.20 

0.75 

0.94 



(0.36) 

(0.33) 


(0.47) 

(0.33) 

(0.37) 


Table 2. Pairwise nucleotide diversity (in per cent, above diagonal) and net nucleotide diversity (in per 
cent, below diagonal) between study units (1-5, see Table 1) of Asian elephant (standard error in pa¬ 
rentheses). 



1 

2 

3 

4 

5 

1 

- 

0.90 (0.37) 

1.18 (0.47) 

0.88 (0.37) 

1.19 (0.45) 

2 

0.06 (0.03) 

- 

1.21 (0.45) 

0.90 (0.40) 

1.06 (0.43) 

3 

0.13 (0.06) 

0.22 (0.09) 

- 

1.17 (0.47) 

1.38 (0.55) 

4 

0.06 (0.03) 

0.18 (0.09) 

0.22 (0.11) 

- 

1.26 (0.53) 

5 

0.28 (0.10) 

0.21 (0.09) 

0.31 (0.12) 

0.41 (0.17) 

- 


Discussion 


Cytochrome b variation in the Asian elephant 

Data available on sequence variation of the mitochondrial cytochrome b gene have 
proved to be inconsistent with respect to genetic variation within mammalian populations 
(see e. g. Shields and Kocher 1991; Randi et al. 1994 for bears; Perry et al. 1995 for 
seals; Carr and Hughes 1993 for deer; Wayne and Jenks 1991; Girman et al. 1993 for 
wild canids). 
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Using restriction enzymes for cutting a 2,450 bp fragment of mtDNA, Georgiadis et 
al. (1994) detected considerable genetic variation within and among populations of the 
African elephant (Loxodonta africana). This result was interpreted to indicate a large 
long-term effective population size and high gene flow between subpopulations. Exhibit¬ 
ing a total of 7 variable sites and considerable haplotype diversity within each of the local 
populations studied, also the Asian elephant, assumed to originate from Africa and to 
have invaded southern and southeastern Asia during the Pliocene (cf. Haynes 1991), can 
be considered at least moderately polymorphic. Despite a remarkable decline and increas¬ 
ing isolation of local Asian elephant populations in historical times due to human exploi¬ 
tation (Kurt et al. 1995), the observed variation may indicate the ancient effective 
population size to have been large. 


Relationships among haplotypes in comparison with their geographic distribution 

Haplotype MAX VI was previously described in six Asian elephants without mentioning 
their geographic origin (Hagelberg et al. 1994). This haplotype is the one most closely re¬ 
lated to those of both the African elephant and the Siberian woolly mammoth, and taking 
also into account its widespread geographic distribution it must be considered a very bas¬ 
ic one. According to their geographic distribution the same holds for haplotype MAX V 
and, although being absent or rare in the central and easternmost part of our study area, 
respectively, for haplotypes MAX III and MAX II. Haplotype MAX I, apparently con¬ 
fined to Sri Lanka may have locally arisen through a new mutation, which is in agreement 
with its marginal position in the median graph. However, due to their more central posi¬ 
tion in the median graph, the other haplotypes found only in single study units (MAX IV, 
MAX VIII) or in a group of adjacent study units (MAX VII) may rather represent relics 
of a formerly more widespread polymorphism. Altogether, the disjunct distribution of 
many haplotypes suggests that they had already been present at the time when elephant 
colonized southern and southeastern Asia during the Pliocene (c. f. Haynes 1991). This ar¬ 
gument is supported by the presence of several well differentiated haplotypes at Cyt b, a 
generally rather conservative gene with respect to variation at the population level 
(cf. Avise 1994). The alternative hypothesis, a segregation of haplotypes MAX IV and 
MAX VIII in the Indian subcontinent, and of haplotype MAX VII in Southeastern Asia, 
respectively, from a basic set consisting of haplotypes MAX VI, MAX V, MAX II, and 
MAX III in the original founder population would imply a considerable extent of homo- 
plasy, which is unlikely to have arisen at the population level. 


Genetic diversity within populations 

The assumption of local losses of Cyt b haplotypes in Asian elephant is supported by lev¬ 
els of nucleotide and haplotype diversity found in extant populations. However, it is inter¬ 
esting to note that, in a comparison of populations studied, estimates of nucleotide 
diversity do not correspond with the respective levels of haplotype diversity. This can be 
explained by both the number of haplotypes and the extent of nucleotide divergence 
among haplotypes being involved in the calculation of nucleotide diversity. The largest 
number of haplotypes is found in Sri Lanka, but with the exception of two animals with 
haplotype MAX VI all of them belong to group A. Hence, the comparatively high nucleo¬ 
tide diversity in that area mainly reflects a large number of haplotypes present. In north¬ 
ern India/northern Myanmar the number of haplotypes found is smaller than in Sri 
Lanka, but the proportions of haplotypes belonging to group A and B, respectively, the 
two being separated by two mutation steps, are more similar than in the previous case. In 
eastern Thailand/Vietnam the high nucleotide diversity is almost exclusively due to a 
large nucleotide divergence among only a few haplotypes. As a methological consequence 
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we suggest that both haplotype diversity and nucleotide diversity are essential for a prop¬ 
er characterization of genetic variability in populations. 

The most striking result with respect to levels of genetic variability is the high haplo¬ 
type diversity present on the island of Sri Lanka as compared to the respective values 
found in most areas of the Indian subcontinent and in southeastern Asia. During the last 
five centuries, nowhere in the range of the Asian elephant has capturing and killing by 
Mogul and colonian forces been carried out as extensively as in the area between south¬ 
ern and northeastern India. There numbers of elephants have declined from at least 
375,000 at 330 BC to less than 20,000 at present (i. e. to about 5 per cent of the original 
population). In southern Myanmar and in Thailand large numbers of elephants had been 
caught for timber extraction between the late 18th and the early 20th century (Evans 
1910; Stracey 1963). In addition it must be stressed that, due to habitat destruction by 
man, on the Indian subcontinent the remaining elephant population is subdivided into a 
number of mostly small isolates, which are exposed to genetic drift (Kurt 1992, Kurt et 
al. 1995). A similar situation prevails in the region from southern Myanmar to Vietnam 
(Santiapillai and Jackson 1990). 

By contrast, in Sri Lanka still some 20 per cent of the original elephant population of 
about 12,000 individuals are present, and human encroachment hardly prevented migra¬ 
tion across the whole island until the middle of this century (Kurt 1992). Also in north¬ 
eastern India and northern Myanmar, harbouring the second highest haplotype diversity, 
human impact on population size of elephant had been comparatively low. Altogether, 
the different levels of genetic variation found in populations of Asian elephant are in 
good agreement with records on population histories. 


Genetic diversity among populations 

Although there seems to be a tendency of haplotypes belonging to group A to decrease 
from the western to the eastern part of our study area, haplotype frequencies and esti¬ 
mates of net nucleotide diversity among samples do not correspond very well with the 
geographic distribution of the local populations studied. A certain bias due to limited 
sample sizes of individuals cannot be excluded, but in general this result seems to reflect 
a considerable extent of genetic drift caused by human impact on population size and 
structure. 


Implications for systematics and conservation of Asian elephant 

Estimates of net nucleotide diversity among populations do not reveal a major separation 
of the Sri Lanka population from all other study units. Thus, our data do not provide sup¬ 
port for subspecies distinction between elephants on Sri Lanka and on the Asian main¬ 
land, which is in accordance with the electrophoretic data of Hartl et al. (1995). In the 
light of our present results the apparently fixed difference at the Sod-2 locus reported by 
Nozawa and Shotake (1990) can be explained by short-term drift processes in local popu¬ 
lations rather than by a long-term divergence between subspecies. 

Regarding breeding of Asian elephants in zoos differentiation among wild-living popu¬ 
lations, especially between those of Sri Lanka and southern India, formerly believed to re¬ 
present different subspecies, seems to be low enough to justify cross-breeding of parents 
of geographically different origin. Future management of wild-living and captive popula¬ 
tions in Asia should follow the traditional Sri Lankan practice, where the off-take from 
the wild-living population for maintaining considerable captive stock has been as high as 
necessary but low enough for preserving the highest level of genetic variability as yet 
known (see also Kurt et al. 1995). 


Population genetics and systematics of Asian elephant 


293 


Acknowledgements 

Collection of samples was made possible by H. E. the Minister of Forestry Lt. Gen. Chit Swe (Yangon); 
Lyn de Alwis, Chairman of the Asian Elephant Specialist Group of the TUCN (Colombo and Singapur); 
Dr. (Mrs.) Jayanthi Alahakoon (Colombo); Dr. Visit Arsaithamkul (Bangkok); Prof. J. V. Cheeran 
(Trichur); Dr. C. Hagenbeck and K. Kock (Hamburg); Dr. S. K. Karunarathe (Colombo); Dr. Along- 
KORN Mahannop (Bangkok); U Myat Win (Yangon); Dr. G. Ruempler (Munster); Prof. K. C. Panicker 
(Trichur); Dr. Wisid Wichasilpa (Bangkok). This research was supported by the R. and R. Schlageter 
Foundation (Zoo Zurich), and by A. Haufellner (European Elephant Group, Munich). 

Zusammenfassung 

Populationsgenetik und Systematik des Asiatischen Elefanten (Elephas maximus): Eine Studie aufder 
Grundlage von Sequenzvariation am Cyt-b-Gen von PCR-amplifizierter Mitochondrien-DNA aus 

Haarwurzeln 

Um die genetische Variabilitat und Differenzierung bei Asiatischen Elefanten {Elephas maximus) zu 
charakterisieren, wurde bei insgesamt 53 Individuen aus Sri Lanka, dem Siiden und Nordosten Indiens, 
dem Norden und Siiden von Myanmar (Burma), dem Norden und Osten von Thailand sowie aus Viet¬ 
nam die Sequenzvariation am mitochondrialen Cytochrom-b-Gen untersucht. Ein ca. 480 bp langer Ab- 
schnitt des Gens wurde mittels PCR aus Haarwurzeln amplifiziert und im Thermocycler direkt 
sequenziert. Nach Auftrennung in einer Direkt-Blot-Sequenzieranlage wurden 335 bp analysiert. Sie- 
ben Sequenzpositionen waren polymorph, und acht Haplotypen wurden definiert, die sich durch eine 
bis sechs Transitionen an der dritten Codon-Position unterschieden. Die Haplotypen wichen im Ver- 
gleich zu publizierten Sequenzen an mindestens 18 Positionen vom Afrikanischen Elefanten (Loxo- 
donta africana) und an mindestens 8 Positionen vom Mammut {Mammuthus primigeniiis) ab. 
Entsprechend ihrer Nukleotiddivergenz lieBen sich die Haplotypen zwei Gruppen zuordnen, die sich in 
ihrer geographischen Verbreitung unterschieden. Die Divergenz zwischen Elefanten von Sri Lanka und 
dem asiatischen Festland, die bisher als eigene Unterarten beschrieben worden waren, war zum Teil 
kleiner als zwischen Elefantenbestanden auf dem asiatischen Festland. Die untersuchten Bestande zeig- 
ten deutliche Unterschiede hinsichtlich der Haplotypen- und Nukleotiddiversitat. Die phylogeogra- 
phischen Muster der Haplotypen legen nahe, daB die Asiatischen Elefanten zur Zeit ihrer pliozanen 
Einwanderung nach Siid- und Sudostasien eine zusammenhangende Population gebildet haben. Unter 
Beriicksichtigung der jeweiligen Bestandesgeschichte laBt sich die gegenwartige geographische Differ¬ 
enzierung als Folge einer unterschiedlich starken menschlichen EinfluBnahme auf die BestandesgroBen 
und -strukturen interpretieren. 
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